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On the gamma-ray spectra radiated by protons accelerated in SNR shocks 
near molecular clouds: The case of SNR RX J1713.7-3946 


M.A. Maikov 1 , P.H. Diamond 1 and R.Z. Sagdeev 2 
Abstract 

Cosmic rays (CRs) are thought to be accelerated in SNRs. The most favorable sit¬ 
uation for proving that the main, hadronic CR component is accelerated there is when 
CRs interact with dense gases, such as molecular clouds (MC) which surround the SN 
shock. Here, a new mechanism of spectrum formation in partially ionized gases near 
SNRs is proposed. Using an analytic model of nonlinear diffusive shock acceleration, 
we calculate the spectra of protons and estimate the resulting y-ray emission occurring 
when the SNR shock approaches a MC. We show that the spectrum develops a break 
in the TeV range and that its GeV component is suppressed. These modifications to 
the standard theory occur because of the proximity of the partially ionized MC-gas and 
because of the physics of particle and Alfven wave propagation inside the gas. Possi¬ 
ble applications of the new spectra to the recent CANGAROO and HESS observations 
of the SNR RX J1713.7-3946 are discussed. 

Subject headings: acceleration of particles—cosmic rays —shock waves—supernova 
remnants—turbulence 


1. Introduction 

Cosmic rays (CRs) travel to us through the chaotic magnetic field of the Galaxy. Therefore, 
the only way to connect them to their accelerator is through the on-site radiation. Accelerated 
electrons have already been detected in supernova remnant (SNR) shocks (Koyama et al. 1995). 
However, electrons comprise only 1-2% of the CR intensity above 4-5 GeV. Therefore, distin¬ 
guishing between leptonic and hadronic origin of the observed TeV emission from SNRs is a key 
to the proof of the supernova origin of CRs. Simultaneous monitoring of X- and y— energy bands 
is an indispensable tool for that purpose. The same TeV electrons radiate in X- and y-rays via syn¬ 
chrotron and inverse Compton (IC) mechanisms, respectively {e.g., Stumer et al. 1997; Reynolds 
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1998). Accelerated protons, can be detected only near their own (most interestingly TeV) energy 
band through their interaction with an ambient gas. Therefore, the “smoking gun” for proton ac¬ 
celeration should be the y-ray emission, without an X-ray emission that can be identified as the 
electron synchrotron radiation. Molecular clouds (MC) adjacent to the SNR will dramatically 
enhance proton visibility (e.g., Aharonian et al. 1994; Drury et al 1994). 

However, to conclusively detect acceleration of super-TeV nuclei turned out to be a very dif¬ 
ficult task. Despite extensive search campaigns (e.g., Buckley et al. 1998; Volk 2000), only three 
SNRs have shown detectable TeV emission so far (Tanimori et al. 1998; Muraishi et al. 2000; 
Aharonian et al. 2001; Enomoto et al. 2002). Nevertheless, a signature of protons accelerated to 
super-TeV energies was reportedly discovered in the SNR RX J1713.7-3946 by the CANGAROO 
team (Enomoto et al. 2002). Recently, this remnant has been confirmed as a TeV source, with 
significantly reduced systematic and statistical errors in the range 1-10 TeV, by the HESS collab¬ 
oration (Aharonian et al. 2004). This remnant is of a shell type, typical for the major acceleration 
models, and rather abundant in the Galaxy. Therefore, after some verification of the data and 
analysis, it could serve as a direct evidence of the CR - SNR connection. 

A number of groups (e.g., Butt et al. 2002; Reimer and Pohl 2002) reanalyzed the data of 
Enomoto et al. (2002) and claimed that the nucleonic interpretation is “highly unlikely” because 
of an alleged inconsistency with the commonly assumed first order Fermi acceleration theory. We 
argue below that if (as it is assumed for this SNR by Butt et al. 2001; Enomoto et al. 2002), some 
of the accelerated particles begin to interact with a partially ionized dense gas such as a MC in the 
northwestern rim of the remnant, the standard acceleration model is not applicable. In particular, 
the nonlinearity of the acceleration process as well as Alfven wave evanescence in the MC are 
essential. Therefore, the analyses by both sides of the controversy are oversimplified and we 
believe that a conclusive case for or against nucleonic origin of the emission has not been made. 
Moreover, a recent analysis by Pannuti et al. (2003) indicates that the standard acceleration model 
applied to TeV electrons does not adequately fit the TeV data, either, unless the X-ray emitting 
filaments have an unreasonably low filling factor of 10 \ Lazendic et al. (2004) have increased it 
to 10 -2 by allowing a smoother spectrum cut-off. 

These analyses indicate that modifications to the accelerated particle spectra are required. 
Namely, a low-energy cut-off above 1 OOGeV and/or a spectral break would support the nucleonic 
scenario (Reimer and Pohl 2002) while models with an exponential cut-off rather than a break 
provide only a bad fit. Such modifications are suggested below. 

The controversy about the TeV observations of the SNR RX J1713.7-3946 is fundamental 
and will clearly impair understanding of future observations of SNRs nearby MCs. Therefore, in 
this Letter we revisit shock acceleration theory and include the following physical phenomena in 
the model: (i) the nonlinearity of the acceleration process, i.e., the modification of the flow by 
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accelerated particles (ii) a position dependent low-energy cut-off of accelerated particles ahead 
of the shock, which is well known to be present in analytic solutions for the shock acceleration 
problem. Note that this cut-off is usually ignored in shock acceleration models since only the 
downstream solution is considered. When the dense gas is present upstream of the shock this cut¬ 
off should be included in the calculation of y-radiation. (iii) impairment of CR confinement in a 
dense gas (MC) due to non-propagation of some Alfven waves. The latter phenomenon produces a 
break (curvature) in the particle and radiation spectra in the TeV energy range and, in combination 
with (i), results in a spectral slope significantly different from the usual predictions of linear theory 
supplemented with the high energy cut-off, which was used in the analyses of Enomoto et al. 
(2002) and Reimer and Pohl (2002), for example. 


2. Low-energy cut-off and the absence of thermal Jt-ray emission 

The diffusive shock acceleration (DSA) mechanism operates by scattering charged particles 
off magnetic irregularities (Alfven waves) across a shock, such as a SN blast wave (e.g., Blandford 
and Eichler 1987). In the linear acceleration regime, the momentum distribution is a power-law, 
P 2 fo{p ) = Cp~ s where C is a normalization constant and the index 5 = (r + 2)/(r — 1) depends 
only on the shock compression ratio r. As soon as the flow is disturbed by the pressure of ac¬ 
celerated particles, the spectrum becomes more complex. It flattens toward higher energies and 
may become steeper at lower energies. The spectrum ff (p) is the downstream spectrum, which 
is coordinate independent in both the linear and nonlinear regimes. Here the upstream spectrum 
is more important since the interaction with the adjacent MCs starts upstream where a low energy 
cut-off occurs. Indeed, the solution upstream, valid for both the linear and nonlinear regimes, reads 
(e.g., Maikov and Drury 2001) 


f(p,x) =/o(p)exp 


- [ X udx 
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Here fo(p) is the spectrum downstream, a(p) ~ 1. The a; coordinate points upstream (x > 0) from 
the shock (x = 0) and u(x) >0 is the speed of the plasma flowing into the shock (u(x) ~ V s i wc k)- 
Now, k(p) grows with p, most likely linearly, K m k Ip. According to eq.(l), there exists a CR 
precursor, which is of the length Iqr — k (p ma x) / V shock ■ Furthermore, high energy particles diffuse 
ahead of low energy particles so the spectrum (1) has an exponential low-energy cut-off at 
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For the Bohm diffusion we have 
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P min/ me « iy shock! c) b m { x ! 10 12 cm) (3) 

where B M is the magnetic field in pG. This low-energy decay can explain why the same MC can 
be visible in the TeV energy range and invisible in the GeV range. Let the near-most edge of an 
adjacent MC be upstream of the shock at a distance x = xmc, such that 1 GeV/c < p m in( x MC) < 
Pmax ■ This low-energy cut-off p m i„ is shown in Fig.l 1 that will be discussed below. The last 
condition means that the leading edge of the CR precursor (filled with TeV protons) has already 
penetrated into the MC, but the sub-shock itself (with the GeV particles ahead of it) still has not. 
Obviously, the CANGAROO and HESS could detect the proton TeV emission while the EGRET 
could not, because of the low density inside the wind bubble the sub-shock is located in. The strong 
X-ray emission expected when the sub-shock crashes into the dense MC will not yet be visible, 
either. 


3. High-energy spectral break 

In the presence of weakly ionized dense gas the particle and emission spectra undergo signif¬ 
icant modifications in the TeV energy range, where a spectral break can form. According to Eq. 
(1) the accelerated particles occupy an extended precursor of the size 

IcR ~ ^{Pmax)/yshock A r g{jPmax)Cshock (4) 

ahead of the shock. Here r g is the particle gyroradius. One can estimate Icr as 
Icr ~ ( Pmax/mc ) {c/Vshock) B u ' 10 12 cm. Thus, the shock precursor may be as long as 10 18 cm for 
Pmax ~ lOTeV. Therefore, the accelerated particles start to interact with a MC before it becomes 
significantly ionized by the shock wave or the ionizing precursor (Drane and McKee 1993). Hence, 
they propagate in the MC under conditions of strong ion-neutral collisional damping of the self¬ 
generated Alfven waves that are needed to confine accelerated particles. More importantly, there is 
a gap in wavenumber space at k\ < k < kj where the waves do not propagate (Kulsrud and Pearce 
1969; Zweibel and Shull 1982). Here, it is critical to realize that waves literally do not exist , as 
opposed to simply being damped by collisions. The above wave evanescence range is bounded by 
k\ — vi n /2 Va and G = 2-y/pj/PoV in/V a, where p,/po 1 is the ratio of the ion to neutral mass 
density, v,„ is the ion-neutral collision frequency and Va — B/sJ 47tp7 is the Alfven speed. The 
resonance condition for the wave generation and particle scattering off them is kp»/m = ±(fl c , 


*We have chosen the following values for the parameters in eq.(3): V s hock = 1000 km/sec, B u = 3, and the distance 
from the sub-shock to the MC, x = xmc ~ 7 • IO'V/h. Note that for the visibility of the TeV and higher energy protons, 
it is sufficient that xmc < 10 17 cm. 
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where oo c = eB /me and /?| is the parallel (to the magnetic field) component of particle momentum. 
Therefore, particles having p\ < p 1 1 < p 2 , where 


pi = 2 V A m(P c /v in , 




(5) 


are not scattered by any waves, so that the confinement of all particles with |p|| | > p\ is dramat¬ 
ically degraded. We assume that the gap in p\ < p < pj is sufficiently broad, so that resonance 
broadening (e.g., Achterberg 1981) does not bridge it. We also assume that there is no significant 
background turbulence at the scales A, > 2 n/k\, that could reduce particle diffusivity to a limit 
which is not significantly higher than the Bohm value. For the parameters used in our calculations 
below, k\ ~ 10~ 14 cm~ 2 . Then, particles with \p\\ | > p\ will escape from the MC upstream of the 
shock along the magnetic field at a speed comparable to c or, at least, their confinement time to the 
shock will be much shorter than those with |pn | < Pi- Only the latter particles will generate 7t° 
mesons and y emission efficiently in the MC. Their pitch angle distribution for p> p\ is, however, 
limited to the interval |/j| < p\ / p. Therefore, the contribution to the y emission of particles with 
momentum p > p\ is reduced by the phase volume filling factor p cr(f = pi/p■ The resulting emis¬ 
sion spectrum is thus one power steeper as compared to the standard calculations (e.g., Berezinsky 
and Ptuskin 1989; Drury et al 1994; Naito and Takahara 1994) based on the isotropic particle distri¬ 
bution. The latter yields the y spectrum which reproduces that of the particles up to about 0.1 cp max 
where it declines. Hence, the energy spectrum of y emission must be the same as that of the parti¬ 
cles, for the photon energy £ y < £/„- ~ 0. \cp\ (particle momenta p < p\), and it must scale 

as £y 5-1 for £y > £/„■ up to about 0.1 cp max (particle momenta p\ < p < p m ax )• The break energy 
Zbr can be estimated using Eq. (5) for p\ by substituting v,„ = no (cV), where o is the cross-section 
of the ion-neutral collisions and V is the collision velocity averaged over a thermal distribution. 
Using an approximation of Kurlsrud and Cesarsky (1971) for (cV), p\ can be estimated as 


Pi/me ~ 16 B 3 T 4 0 A n 0 1 n i 1/2 . (6) 

where T 4 is measured in the units of \QrK, no and n, (number densities corresponding to the mass 
densities po and p,) -in cm~ 3 . 

To illustrate these results we consider three different acceleration regimes, corresponding to 
the three substantially different shock compression ratios r, shown by the points 1-3 in the inset to 
Fig.l. (Here r means a total shock compression, that is adiabatic compression across the precursor 
times the subshock jump.) As it can be seen from this plot, the shock compression is difficult to 
calculate accurately due to its sharp dependence on the injection rate (usually not well known) and 
other parameters, such as the maximum momentum. We discuss this spectrum in more detail in 
the next section. 



4. Discussion and Conclusions 


In this paper we have considered the upstream particle spectrum since this is necessary for 
interpretation of cases where a shock is expanding into a low-density pre-supernova wind bubble 
and is approaching a denser material such as a swept-up shell or a MC. Perhaps the most dis¬ 
cussed SNR of this kind is the RX J1713.7-3946 claimed by the CANGAROO team to be a long 
anticipated proton super-TeV accelerator (Enomoto et al. 2002). Significant part of the TeV emis¬ 
sion detected by the CANGAROO and HESS (Aharonian et al. 2004) instruments comes from the 
northwestern rim of the remnant where the interaction with a molecular cloud is believed to take 
place (e.g., Slane et al. 1999; Hiraga et al. 2004). Therefore these observations are pertinent to the 
subject of this letter. 

To demonstrate that the observed TeV spectra are consistent with the mechanisms suggested in 
this paper, we put the calculated upstream proton spectra on one plot with the CANGAROO/HESS 
and EGRET photon fluxes that would be radiated by these protons. Since the emission generated 
by protons has typically an order of magnitude higher energy, we up-shifted the observed emission 
energy accordingly in Fig. 1. As the spectrum between GeV/c < p < p\ is considerably flatter 
than in the case of the linear theory, its magnitude at p — p\ is an order of magnitude higher for 
the same injection rate. This lowers the ambient gas density required to produce the same TeV 
emission as computed by Enomoto et al. (2002) using the linear spectrum. In addition, this allows 
to accommodate the EGRET points, should they belong to the same object (Reimer and Pohl 2002). 
However, because of the controversy about the distance to the remnant (6 vs 1 kpc) (Slane et al. 
1999; Fukui et al. 2003) and other uncertainties, it is impossible to present detailed normalized fits 
of the photon flux in this short letter. 

The form of the spectrum is marked by increasing the spectral index by one at the break 
momentum p = p i, Fig. 1. This is crucial to fitting the observed spectra without imposing a 
specific form and value of the energy cut-off. The latter procedure would give only a poor fit 
(e.g., Reimer and Pohl 2002). Note that there are two CANGAROO points that are not in a good 
agreement with both the HESS and with the theoretical spectrum. It should be noted, however, 
that any disagreement between the three sets of spectral points is substantially enhanced in this 
particular spectrum format (the particle phase density is multiplied by p 1 ). Except for these two 
points the agreement is remarkably good. Interestingly enough the agreement is excellent for a 
subset of HESS points with the lowest and highest energy points excluded. 

Unlike the power-law index, it is more difficult to constrain the parameters determining the 
position of the break on the spectrum given by Eq. (6) due to the poor information about the 
target gas for p — p reactions (MC). For the particular case shown in Fig.l, we have chosen such 
combination of parameters in Eq. (6) that p\ amounts to p\ ~ l.STeV/c. MCs in general are 
known to be “clumpy” with the inter-clump gas density of 5-25 cm -3 and a less than 10% filling 
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factor (see, e.g., Chevalier 1999; Bykov et al. 2000, and references therein). The results shown in 
Fig. 1 are obtained for the following values of parameters in eq.(6): B^ « 14, T 4 « 0.01, no & 23, 
and ni ~ 0.23 which yields pi/mc ~ 1.8 ■ 10 3 . By contrast, in the model of Ellison et al. (2001) 
much lower target density is assumed since its authors do not consider the possibility that the CR 
precursor reach the dense gas while the shock itself is still in the wind bubble. Therefore they 
concluded that the p — p reactions do not contribute significantly to the TeV emission because 
of the lack of target protons. In contrast to the conclusion of Ellison et al. (2001) and to that of 
Enomoto et al. (2002) the most recent HESS observations allow their authors (Aharonian et al. 
2004) to assume that both protons and electrons contribute to the TeV emission significantly. 

Note that a complementary mechanism of suppression of the low-energy (GeV) emission has 
been suggested earlier by Aharonian and Atoyan (1996) and discussed recently in the context of 
SNR RX J1713.7-3946 by Uchiyama et al (2003). This mechanism requires an impulsive release 
of accelerated particles at some distance from a MC (target) and subsequent (energy dependent) 
diffusive propagation of CRs to the target. By contrast, the mechanism discussed in this paper is 
based on a quasi-stationary solution of the acceleration problem (Eq. [1]) which implies that the 
accelerated particles are bound to the propagating shock front via self-generated Alfven waves. It 
is clear, however, that in both cases the suppression of the low-energy emission is based on a slower 
diffusion of low-energy particles. Note that Lazendic et al. (2004) also discussed qualitatively the 
same effect. It is also important to mention that GeV emission may not be suppressed all over 
the remnant, and was demonstrated by Butt et al. (2001) to be likely of proton origin. Obviously, 
the spectral modifications considered in this letter are equally applicable to electrons of similar 
rigidity. 

It should be pointed out that our model in its current version does not provide a complete 
fit to the data, particularly in the X- and radio energy bands. The TeV-GeV fluxes are given in 
arbitrary units, so only the form of the spectrum (including possible GeV signal) is provided here. 
A complete fit would require the convolution of the proton spectrum with appropriate emissivities 
given the densities and the distance. The thermal peak in Fig.l is a Maxwellian, approximately 
normalized relative to the high-energy part of the spectrum, in accordance with the calculated 
downstream temperature without turbulent heating. More accurate calculation of the thermal peak 
would require adequate injection model and turbulent gas heating in the precursor {e.g., Berezhko 
et al. 1996; Maikov and Drury 2001; Kang et al. 2002). Such calculations of the thermal X-ray 
emission should be tested specifically against very low observed upper limits. One potentially 
observable prediction of our model is the break in the pitch angle averaged particle spectrum at 
p — Pi that would be produced by almost instant losses along the field lines of all particles with 
|p|| | > pi once they enter a partially ionized gas with significant ion-neutral collisions. This break 
may very well be in the currently “obscure” energy range between the EGRET and the ground 
based Cherenkov telescope energy bands. The next generation of space gamma telescopes, such 
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as GLAST, will be able to explore this energy range. 

The principal results of this letter, which are the spectrum softening by one degree above 
the spectral break, and the possibility of a low energy cut-off clearly show that these physical 
phenomena need to be included in the models to conclusively differentiate between the nucleonic 
and leptonic sources of the TeV emission from the remnants nearby molecular clouds such as the 
RX J1713.7-3946. Previous claims to the contrary have ignored necessary refinements in the DSA 
theory beyond the level of test particle theory. 


We thank Felix Aharonian and Heinz Volk for helpful discussions as well as an anonymous 
referee provided useful comments on the manuscript. Support under NASA grant ATP03-0059- 
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Fig. 1.— Inset: response of the shock structure to the particle acceleration. The flow compression 
ratio r is shown as a function of particle density uqr injected into the acceleration from the ther¬ 
mal distribution calculated for two maximum momenta p max and the shock Mach number M=80. 
The injection parameter v ~ (cp m j/mV^ ho/ . k )ncR /« 1 , where n \ is the plasma density upstream of 
the shock. The details of these calculations and those of the spectra shown in the main figure can 
be found in (Maikov and Drury 2001). Main figure: spectra of accelerated particles behind the 
shock front, for p max /mc = 10 5 . Three solutions correspond to the points 1-3 on the compression- 
injection diagram in the inset. Note that the test particle (linear) spectrum would be a horizontal 
line. The break energy (momentum) is at p/mc — 1.8 • 10 3 . The CANGAROO-II (Enomoto et al. 
2002) and HESS points are adopted from Aharonian et al. (2004) while EGRET (Hartman et al. 
1999) points- from Reimer and Pohl (2002). Vertical lines indicate: the injection momentum (sep¬ 
arating thermal and nonthermal particles), the 1.2 GeV proton energy threshold of 7t° production, 
the possible low-energy cut-off p m in(x ) and the spectral break p\ (see text). 







